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Abstract. We consider a general class of Glauber dynamics reversible with respect 
to the standard Ising model in 1 d with zero external field and inverse temperature (3 
strictly larger than the critical value f3 c in dimension 2 or the so called "slab threshold" 
j3 c in dimension d ^ 3. We first prove that the inverse spectral gap in a large cube 
of side N with plus boundary conditions is, apart from logarithmic corrections, larger 
than N in d — 2 while the logarithmic Sobolev constant is instead larger than iV 2 in any 
dimension. Such a result substantially improves over all the previous existing bounds 
and agrees with a similar computations obtained in the framework of a one dimensional 
toy model based on mean curvature motion. The proof, based on a suggestion made by 
H.T. Yau some years ago, explicitly constructs a subtle test function which forces a large 
droplet of the minus phase inside the plus phase. The relevant bounds for general d > 2 
are then obtained via a careful use of the recent L 1 -approach to the Wulff construction. 
Finally we prove that in d — 2 the probability that two independent initial configurations, 
distributed according to the infinite volume plus phase and evolving under any coupling, 
agree at the origin at time t is bounded from below by a stretched exponential exp(— Vi), 
again apart from logarithmic corrections. Such a result should be considered as a first 
step toward a rigorous proof that, as conjectured by Fisher and Huse some years ago, 
the equilibrium time auto-correlation of the spin at the origin decays as a stretched 
exponential in d = 2. 
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1. Introduction 

In a finite domain, the reversible Glauber dynamics associated to the Ising model relaxes 
exponentially fast to its equilibrium measure. Nevertheless, this simple statement hides 
a wide range of behaviors depending on the temperature, the domain and the boundary 
conditions. 

In the uniqueness regime (when the temperature is large enough), the speed of relax- 
ation is uniform with respect to the domains and the boundary conditions. We refer to 



Martinelli [Ma] for a complete account of this theory. The occurrence of phase transition 
drastically modifies the behavior of the dynamics and new physical features slow down 
the relaxation; among those, the nucleation and the interface motions. Metastability is 
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characteristic of these slow phenomena since the system is trapped for a very long period 
of time in a local equilibrium. In this case, the relaxation mechanism is so slow that 
the time of nucleation can be expressed in terms of equilibrium quantities. In particular, 
it was proven by Martinelli (see eg. |Ma| and references therin) that for free boundary 
conditions the asymptotic of the spectral gap with respect to the size of the domains is 
related to the surface tension and the main mechanism driving the system to equilibrium 
is nucleation of one phase inside the other. A complete picture of the nucleation process 
in 1? in the framework of metastability was obtained by Schonmann and Shlosman in [3S2]. 



In this paper, we are interested in a different regime in which the relaxation to equilib- 
rium is driven by the slow motion of the interfaces. This is the case of the Ising model 
in a large box with plus boundary conditions. When a droplet of the minus phase is 
surrounded by the plus phase, it tends to shrink according to its curvature under the 
action of the non-conservative dynamics on the spins close to the interface. This subtle 
phenomenon has been studied rigorously only in rare instances: by Spohn |f5p|| in the 



case of Ising model at zero temperature (see also Rezakhanlou, Spohn [RS]), by Chayes, 



Schonmann, Swindle [CSS] for a variant of this model and by De Masi, Orlandi, Presutti, 



Triolo [DOPT1, DOPT2 for the Kac- Ising model. Notice also that the motion by mean 



curvature plays a key role in the coarsening phenomenon, as it has been shown recently 



by Fontes, Schonmann, Sidoravicius [FSS|. For positive temperatures, a mathematical 



derivation of similar results seems to be more challenging. 

A way to capture some insights into the slow relaxation driven by interface motion is to 
estimate spectral quantities related to the generator of the Glauber dynamics. We prove 
that for any dimension d ^ 2, in the phase transition regime and with plus boundary 
conditions, the logarithmic-Sobolev constant for a domain of linear size N diverge at least 
like iV 2 (up to some logarithmic corrections). This can be considered as a first characteri- 
zation of the slow down of the dynamics and is in agreement with the heuristics predicted 
by the motion by mean curvature. In the same setting but d = 2, we prove that the inverse 
of the spectral gap grows at least like N (up to logarithmic corrections). In dimension 
d ^ 3 our argument fails to produce a result on the divergence of the spectral gap. 

Let us stress that we have not been able to derive matching upper bounds; the best 
existing bounds have been proved only in d = 2 and are of the form exp(\/iV (log iV) 2 ) (see 
Y W| ] ) . However, an exact computation for a toy model based on mean curvature motion 



seems to confirm that the polynomial asymptotics we obtain are correct (see Section [/]). 
The proof boils down to bound the variational formula for the Poincare and the Log- 
Sobolev inequalities by choosing an appropriate test function. This reduces the problem 
to a computation under the equilibrium Gibbs measure. The main difficulty is to recover 
polynomial bounds by using only the exponential estimates provided by the equilibrium 
theory of phase segregation (see |BIV | and references therein). This is achieved by the 



choice of a subtle test function which was suggested some years ago by H.T. Yau. 

The second part of the paper (section ||) applies the result on the lower bound on the 
inverse of the spectral gap to investigate the relaxation in the infinite domain Z 2 . Thanks 
to an heuristic argument based on the motion by mean curvature, Fisher and Huse flHFH 
conjectured that the equilibrium time auto-correlation of the spin at the origin decays as 
a stretched exponential exp(— \/t) in d = 2. We provide a first step towards a rigorous 
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proof by showing that a dynamical quantity strictly related to the auto-correlation cannot 
relax faster than exp(— Vt). 



2. The model and the Main Results 

In this section we define the model and fix some useful the notation, recall some basic 
facts about the Ising model below the critical point and finally state our two main results. 

2.1. The standard Ising model. Let A be a generic finite subset of Z d , with d ^ 2. 
Each site i in A indexes a spin crj which takes values ±1. The spin configurations {<7j}jgA 
have a statistical weight determined by the Hamiltonian 



2 

|i-j|=l l»-j|=l 



over 



i u 3 i 



where a = {(Ji}i e \c are boundary conditions outside A. 

The Gibbs measure associated to the spin system with boundary conditions a is 

= {a x } xeA , /4(<x) = -L- exp (-/?#>)) , 

P,A 

where (3 is the inverse of the temperature (/? = y) and is the partition function. If 
the boundary conditions are uniformly equal to 1 (resp. —1), the Gibbs measure will be 
denoted by p\ (resp. /i A ). 

The phase transition regime occurs at low temperature and is characterized by sponta- 
neous magnetization in the thermodynamic limit. There is a critical value (3 C such that 

V/? > 0c, lim n+(a ) = - lim /x7(<r ) = m* > . (2.1) 

A^Z d A^Z d 

Furthermore, in the thermodynamic limit the measures //^ an d A*a converge (weakly) to 
two distinct Gibbs measures fi + and [i~ which are measures on the space {±l} zd . Each of 
these measures represents a pure state. In dimension d ^ 3, we also denote by (3 C ^ 0, 



■c 



the "slab critical point" (see [ACCER] and |P|]) which is conjectured to coincide with (3 C . 
For convenience we set (3 C = (3 C in dimension 2. Our proofs rely on results of equilibrium 
phase coexistence for the Ising model which are restricted to values (3 > (3 C (for technical 
reasons) . 

The next step is to quantify the coexistence of the two pure states defined above. Due 
to the lattice structure, the surface tension is an-isotropic. Let A = {— N, . . . , N} d , let 
n be a vector in § rf_1 such that n ■ e \ > and let a be the following mixed boundary 
conditions 

\/i £ A c , o; = 



if ft - i ^ 0, 
if n ■ i < 0. 



The partition function with mixed boundary conditions is denoted by Z^ N (n) and the 
one with boundary conditions uniformly equal to +1 by N . 
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Definition 2.1. The surface tension in the direction n 6 §* d 1 , with n ■ e\ > 0, is defined 
by 

r(n)= lim -^log^^. (2.2) 



We refer to Messager, Miracle-Sole and Ruiz [MMR] for a derivation of the thermody- 
namic limit (|2.2| ), Associated in a natural way to the surface tension is the Wulff shape 
which describes the optimal shape of a droplet of the minus phase immersed in the plus 
phase. 

Definition 2.2. The Wulff shape is the convex set in M. d given by 

W= P| {x€R d ; x ■ n ^ r(n)\ . (2.3) 

rigs'*- 1 

The Wulff shape with volume 1 is denoted by W d . Finally in what follows we will choose 
for simplicity the finite set A as the domain Wjy = A r W d nZ rf , instead of a cube of side N. 
The corresponding Gibbs measure on Wtv with + boundary conditions will be denoted by 
Mat- 



2.2. The Glauber dynamics. The stochastic dynamics we want to study is defined by 
the Markov generator given by 

(£+/){*) = £ 4(a)V x f(a) 

zeWiv 

where the values of a outside Wat are kept fixed identical to +1 and S/ X f(a) = [f(cr x ) — /(<r)]. 
On the flip rates c x (a) we assume 

(i) k~ l < c x (a) < k for some k and any x, a 

(ii) reversibility w.r.t. the Gibbs measure fit, 

(iii) finite range 

Remark 2.1. It is possible to check (see e.g. |L| or |Ma| j that it is possible to extend the 
above definition of the generator directly to the whole lattice Z* d and get a well defined 
Markov process on Q := {0, l} z<i . We will refer to the latter as the infinite volume Glauber 
dynamics. 

The Dirichlet form associated to takes the form 

4(/,/)= E /4(^)|v,/| 2 ) 

xew N 

and, thanks to assumption (i) on the flip rates it is uniformly bounded from above and 
from below by 

ti( E |V,/| 2 ) :=/,+ ( |V/| 2 ) 

Two key quantities measure the time scale on which relaxation to equilibrium occurs. The 
first one, denoted by SV, is the inverse of the spectral gap of the generator, while the other 
one the logarithmic Sobolev constant L^. They are both characterized by a variational 
principle in that they are the optimal constants in the Poincare inequality 

»Mf,f)<c£+(f,f), v/aW 
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and in the logarithmic Sobolev inequality 

^+(/ 2 log f)<c£+(f ',/), V/GL 2 ^) with^+(/ 2 ) = l 
respectively. As it is well known the quantity Sn measures the relaxation time in an 
' 4) ' 



L 2 (dfi~jy) sense while Ln measures the relaxation time in an L°° sense (worst case for the 



initial condition). More precisely, if pj- + ' N ^ denotes the Markov semigroup generated by 
Cpf and / is an arbitrary function with Hj^(f) = then 

^([^ (+,iV) /] 2 )<^(/ 2 )exp(-^). 

In many cases e.g. at high temperature the two quantities are of the same order but it 
may very well happen that they are quite different. We will argue later on that the Ising 
model below the critical temperature is actually one of these cases. 

2.3. Main Results. We are finally in a position to state our main results. 

Theorem 2.1. Assume d = 2 and j3 > f3 c . There exists a constant k defending on j3 
such that 

l im (log AT g = +QO r 2A) 

Remark 2.2. As we already pointed out in the introduction, in dimension greater than 
two our choice of the test function to be inserted in the Poincare inequality does not provide 
any non trivial information. 

The next result concerns the large N behavior of the logarithmic Sobolev constant. 

Theorem 2.2. Assume d ^ 2 and (3 > (3 C . There exists a constant k depending on (5 
and d such that 

(logAT Ln = +oo (25) 

N— >oo 1S Z 

Finally we investigate in d = 2 the relaxation in the plus phase for the infinite volume 
dynamics. For this purpose, let us consider an arbitrary coupling of the Glauber dynamics 
in the infinite volume Z 2 . The two processes at time t are denoted by (cr J) (t) , a w (t)) , 
where (r/, lo) are the initial spin configurations. The joint expectation of the process is 
denoted by E. The initial conditions will in general be chosen w.r.t. the product measure 
dfi + (r],Lo) = dn + (r])dfi + (uj), where /i + is the Gibbs measure in the + pure phase. 

Theorem 2.3. There exist positive constants C\,C2 and k independent of the choice of 
the coupling such that 

Vi>0, J dfL + { V ^)t{al{t) ^d%(t)) > C ie xp(-C 2 Vi(logt) K ) . (2.6) 

Remark 2.3. Although we believe that the quantity considered in the theorem is a good 
measure of the time auto-correlation in the plus phase of the spin at the origin, the latter 
is unfortunately only bounded from above by the LHS of fl2.6|) . We have in fact 

^+ ((P t (a ) - m*) 2 ) = fi+ ({Pt(ao) - J d/i + P t (a )) 2 ) 

= / d^( V )((j dv + (u)^(t)-a%(t))Y] ^4 J dfi + ( V ,u;)E{^(t)^^(t)) 
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Remark 2.4. A related result at (3 = +oo was proved recently in [FSS] for the zero 
temperature dynamics (see theorem 1.2 there). 

3. Large Deviations 

In this section we recall some results on the large deviations for the Gibbs measure p^ 
when j3 > f3 c . Our proofs rely on a weak description of phase segregation in terms of 



L 1 -norm. In dimension 2, more precise results can be found in Ioffe, Schonmann [[Sc| and 
in Pfister, Velenik |PVf (for L 1 -concentration statements). The reader is referred to [BTV] 
for a survey on phase coexistence and a complete list of references. 

We consider our microscopic Ising model embedded in W d . LetW N = ^% d n W d and let 
K be a mesoscopic scale (eventually depending on N). The domain W d is partitioned into 
boxes Btv,a", each of them containing K d sites of Wjv: 



j G Z d , Xj = j— £ Wtv, »N,K(Xj) = Xj + 



K K ld 
W 2N 



Let Mx(Nxj) be the microscopic counterpart of Bjv,if (a^), i.e. the sites of inBjv,^^)- 
These boxes are centered on the sites of Wjv,k = fZ d nW t! . As the domain is not regular 
some boxes may not fit inside Wjv, therefore at the boundary we consider a relaxed notion 
of boxes. 

Finally, the local magnetization is defined as a piece- wise constant function on the partition 

Vy € M NtK (xj), M N , K {y) = 1 ^- 1 V Oi. (3.7) 

iff — » 

ieBK(JVaSj) 

The local order parameter Mn,k{u) characterizes the local equilibrium of the mesoscopic 
box containing y. The key result concerning the local order parameters is a trivial con- 
sequence of the results obtained by Pisztora jPj and it is based on the following coarse 
grained description. To each box Bjv^ we associate a mesoscopic phase label UN,K(%j) 
taking values in {—1, 0, 1} 

UN,K(xj) = l{\M N> K(xj)-m*\ ^ \m"} ~ l {\M N , K (xj)+m*\ < \m*} • 

The distribution of the variables {uN : K(xj)} under is dominated by Bernoulli Perco- 
lation. 



Theorem 3.1 ([jPj). For any (5 > (3 C there exists eg > and 7 G]0, 1[ such that the 
following holds uniformly in N: 

V{xi, ...,X£}e W N:K , {un,k(xi) = 0, . . . , u N:K (x e ) = 0) ^ (pk) , (3.8) 
with pk = exp(— cpK 1 ). 



Remark 3.1. For the next results (Propositions |ffJ| , \3.2\ 3.S ) to hold true the mesoscopic 
scale K has to be chosen just large enough ( depending on (3 and on some extra parameter 
5). However in the next sections it will be essential to relate K with the basic scale N 
via the scaling relation K ~ (logiV) 1 / 7 and therefore we will adopt this choice right away 
and denote the corresponding mesoscopic phase labels simply by un- Moreover, since the 
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blocks with label will play an important role in the proof of the main results, they will be 
referred to as the bad blocks. 

In order to state the other results on the large deviations of /i^- we need to introduce 
some more notation. For any 6 > 0, the ^-neighborhood of v G L 1 (W d ) is defined by 

V(v, 5) = {v' G L x (W d ) | \\v' - v\\i < 5} . 

Let O be an open set containing W d . The set of functions of bounded variation in O 
taking values in {—1, 1} and uniformly equal to 1 outside W d is denoted by BV(W rf , {±1}) 
(see [EG] for a review). For a given a > 0, the set of functions in BV(W rf ,{±l}) with 
perimeter smaller than a is denoted by C a . Finally we define the Wulff functional Wp on 
BV(W d , {±1}) as follows. For any v G BV(W d , {±1}), there exists a generalized notion 
of the boundary of the set {v = —1} called reduced boundary and denoted by d*v. If 
{v = —1} is a regular set, then d*v coincides with the usual boundary dv. Then one 
defines 



Wpiy) := / T(n x )dH x , 

Jd*v 

where TC X is the d — 1 Hausdorff measure. The Wulff functional Wp can be extended on 
L^W*) by setting 

W B ( V ) = [ Jrv-rte)***' if v£BV(W d ,{±l}), M 
\ oo , otherwise. 

For any m in [—m*,m*[, the Wulff variational problem can then be stated as, 

mini Wpiy) v G BV(W rf ,{±l}), | / m*v r dr\^m\. (3.10) 
I Jw d J 

If we denote by T> m the set of minimizers of ( 3.10 ) it has been proven by [Ta] that in 
M. d the minimizer is unique up to translations and given by suitable dilation of the Wulff 
shape (2.3). In particular the interfacial energy of W d is given by 

T *=Wp(W d )= f T(n x )dH x . (3.11) 

JdW d 

All that being said the results we are going to use can be summarized as follows. 
Proposition 3.1. There exists a constant C{(3) > such that for any 5 > 

Va > 0, limsup _! log/^ ( U N £V{C a ,5)) < -C(f3)a, 
where V(C a ,5) is the 5 -neighborhood of C a inL 1 (W d ). 

This proposition tells us that only the configurations close to the compact set C a have a 

contribution which is of a surface order. 

The precise asymptotic related to surface tension are 

Proposition 3.2. Uniformly over 5 > 



!im|nf log^i+ (\\u N + < 5) ^ 
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Proposition 3.3. For all v in BV(W , {±1}) such that Wp(v) is finite and for 5 > 

limsu P T7ZiT lo gMjv(ll u iV -^||i < 8) < -Wp{v) +s(6), 
where e(S) vanishes as 5 goes to 0. 

4. The test function 

In this section we define the test function that, if plugged into the Poincare and loga- 
rithmic Sobolev inequalities, will prove theorems ( |2.1| ) and fl2.2| ). As we mentioned in the 
introduction, the form of function described below was suggested by H.T. Yau. 

Fix A G]|t*,t*[, where t* is defined in (3.11). Let g be a smooth non increasing 
function such that 




if s < - 
if s > 



2 



-m 
4 



The mesoscopic scale K is chosen equal to (b log N) t , where 7 was introduced in ([^) and 
6 is a constant which will be fixed later. The test function / has the following expression: 

(\K d \ 
— YsgiM^Kixj))] , VaG{-l,l} w ". (4.12) 

The factor K d stands for the volume of the boxes Mk which equals to K d (with the 
exception of some boxes along the boundary). Notice that / is a non increasing function 
of the spins. 

There are three main features of / that make it quite effective. These are: 

i) The variance of / almost coincides with /x^(/ 2 ) or, put it in another way, /x^(/ 2 ) 3> 

/4(/) 2 ; 

ii) The entropy of f 2 w.r.t. to /x^ is of order N d l ; 

_|_ j 

iii) Let us denote by fj,~L'f the weighted measure d ^ N + = —^rj f 2 where := /^(/ 2 ). 
Then under /Ujy the typical number of non zero terms in |V/| 2 is of the order of 



N' 



1-1 



It is clear that once these properties are established then the proof of theorems 2.1 and 



2.2| should follow quite easily. 

Intuitively the proof of i) , ii) and iii) is based on the following simple heuristic. The 
function / assigns an exponential weight to the configurations with a large number of 
mesoscopic boxes with label un = — 1 because of the choice of the function g. According 
to the large deviation theory, among the configurations favored by /, those with the largest 
ii^r weight form a Wulff droplet of a certain size. Therefore, to compute fJ^(f) or /z^-(/ 2 ), 
we will need to compare, for a given Wulff droplet, the gain given by the exponential factor 
in / or f 2 with the fj,^ probability of creating the droplet itself. It turns out, due to the 
precise choice of the parameter A, that the balance for / is negative and no Wulff droplet 
will appear, while the balance is positive for f 2 and the typical spin configurations under 
fipf will consist of a Wulff droplet of the minus phase of volume ~ N d . That accounts 
for i) and ii). Given the above picture, it is also clear that iii) holds simply because the 
non zero terms in |V/| 2 come only from the bad boxes, again because of the choice of the 
function g. The boundary of the Wulff droplet produces 0(N d ~ l ) of such boxes while the 
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inside of the droplet typically does not contain any bad box because of the choice of the 
mesoscopic scale K. Were K be large but independent of iV then we would always have 
a density of bad boxes and the whole construction would break down. 



4.1. The Variance of /. We are first going to check that 



lim 



lim /4(/ 2 )-/4(/) 2 ^, 
/4(/ 2 ) 



(4.13) 



The function uniformly equal to — 1 in W d is denoted by —1. Let e > 0, then 



where we used the fact that if un(x) = — 1 then Mn,k(x) < — ^n- Proposition implies 
that for ./V large enough 



V+(f) > exp(N d - 1 {2\(l-e)-T*-o(e))) , 



where o(e) vanishes as e goes to 0. 
Next we examine (Hf(f) and prove that 

^ooiV^ 1Og ^ (/) = - 



(4.14) 



(4.15) 



The derivation of an upper bound for requires some technicalities. First fix a 

constant a > where C((3) appears in Proposition 3.1. Then Proposition 3.1 implies 
that for any 5 > 



1 



lim SU P TFd^i lo § ^% (/ 1 {u N ^v(C a ,S)}) < A ~ C (P) a < • 



(4.16) 



jv^oo N d - 

Then fix e > and recall that C a is compact for the L 1 topology. According to Proposition 



3.3 , for 5 small enough the set V(C a , 5) can be covered by a finite union Lr| =1 V(i;j, £$) such 

(4.17) 



that for any i ^ i and N large enough 
1 



— g=iIog^(ujv G V(Vi,Ei)) < -W^(Ui)+e, 

where e$ G (0, e). Noticing that 



IWgVK,^)}) + Mat (/ HuNgViCaJ)}) ) 

and combining ( |4.16| ) with ( |4. 17| ) , we get 

l 

vUf) < E ex P (^"^N - W^(«i) + e(1 + A))) + exp (iV^A - C(P)a 



(4.18) 



i=i 



where denotes the volume of the set {vi = —1} . To check that the spin configuration in 
{un G V(vi,£i)} have a number of blocks of the order of N d \vi\, it is enough to regularize 
Vi by a polyhedral set (see Thm 2.5.1 in | BIV |). 

By the very definition of the Wulff variational problem, for any v G BV(W rf , {±1}) 



(4.19) 
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where we have used the fact that \v\ ^ \W d \ = 1. As A < r*, 

//+(/) ^expjV-^l + A)) +1. (4.20) 



Since e was arbitrary, this implies (4.15). 

Combining (4.14), (|4.15| ) and the fact that 2A > r*, we finally derive (|4.13| ) by choosing e 
small enough. 

4.2. The Entropy of f 2 . We will prove that 

CN d - 1 < ^(log/ 2 ) - lo g/ u+(/ 2 ) , (4.21) 

for a suitable constant C. 

Using the previous strategy, we check that for e > and for TV large enough 



V + N (f) < E ex P (^ 1 (2A|^| - Wp(vi) + e(l + A))) + 1 . 
i=i 

Inequality fl4,19| ) implies 

t 

»nU 2 ) < ^exp(iV d - 1 ((2A-T*)|^|+e(l + A)))+l 



■i=i 

< £exp(iV d - 1 (2A-T* + e(l + A))) +1. (4.22) 
It remains to check that for e > and iV large enough 

4' / ( 1 °g/ 2 ) = ^^' / ( E 9(M N , K (x))) > (l-e)2\N d - 1 . (4.23) 

This is a consequence of the following estimate. For any e > and N large enough 

(1 - o( e ))/4 (/ 2 ) < M + (/ 2 l {uive v(-i, £ )}) • (4.24) 
Let J 1 * = (V(— I, e)) C . First notice that 

sup{-W(u) + 2AM} < sup{M(-r* + 2A)} < (-t* + 2A)(1 - e) . 

We proceed as before and cover the set J-n V(C a , 5) with a finite number of neighborhoods. 
This implies that for any 5 > and N large enough 

Mjv (/ 2 l{ Ujve ^}) < ^exp ((-r* + 2A + (<5))(1 - e)^ 1 ) + 1 . 
On the other hand, 



4(/ 2 l KeV (-i lE )}) > exp((-r* + 2A)iV 



d-l> 



Thus, for N large enough, we derive ( |4.24j ). This implies that 

K d ^(f 2 ( E 5(-^iV,x(x)))) >4(/ 2 (l W£VHe)} ))(l- £ )JV d . 



The inequality fl4.24|) leads to 

K d ^(f 2 ( E 5(-M^W))) ^ /4(/ 2 )(l-e)(l-o(e))iV- 



xGWjv,K 
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Therefore ( pL23|) is complete. Combining ( |Og ) and ( p~23| ), we obtain for A" large enough 

^ J (log/ 2 ) - log/^+(/ 2 ) > ((1 - o(e))2A - (-r* + 2A + + A)))^" 1 - log*, 

^ (t* - o(e)2A - e'(l + A)) A^" 1 - log* . 

For any (3 > (3 C (the true critical point) it is known that t* > 0. Thus, by choosing e and 
e' small enough, we derive ( 4.21| ) for N large enough. 

4.3. The Dirichlet Form of /. The Dirichlet form associated to / can be bounded as 
follows. There is C\ > such that for N large enough 

/4(|V/| 2 ) < C 1 X 2 N d ~ 3 K d n+ (f) = dA 2 iV d - 3 (61ogiV)^+ (/ 2 ) . (4.25) 
By Taylor expansion 

|2 



|V/|' 



E 

ieWjv 



E E iv«/i s 

A 2 



< / 2 (2^||,'||L) E 1 



xeWiv . 



where Qn denotes the number of blocks in Wat with averaged magnetization in [— — 
Using the notation z^'^ introduced in ii) above we can write 



/4(|v/| 2 ) < C1 



A 2 if a 



I'/IlL/^/ 2 )^ 7 ^). 



The estimate (4.25) will follow from the fact that for N large enough 

^ f {Q N ) ^2iV*-i. 

This boils down to check that 

4 J (Qn Hq n >n*-i}) < iV (i exp(-ciV d - 1 ) , 
where c is a positive constant. As p^(f 2 ) 1, we see that 



4 J 1 {Qjv >^-i } J < e X p(2AiV' 1 - 1 ) / u+ (Q* > A r<I_1 J ^2 



rd-1 



(4.26) 
(4.27) 

(4.28) 



Remember that the occurrence of bad blocks is dominated by Bernoulli percolation with 
parameter px = N~ c i lb . Therefore, for b large enough, it is then quite simple to check 
that there is c > such that for A" large enough, 



l4r ( Qn > A^ 1 ) < exp 



cN' 



Combining the previous bound with (4.28), we derive ( 4.27| ). 

Remark 4.1. It would be possible to derive sharper estimates for ( 4.26 ). One expects 

Nevertheless this would not be enough to derive an asymptotic for the spectral gap and 
the Log-Sobolev constant without a logarithmic correction : on finite mesoscopic scales, we 
cannot control the test function. 
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5. Proof of Theorems 2.1 



AND 



2.2 



We are in position to prove the first two main results. 



5.1. Proof of Theorem 2.1. By definition 

When <j) is equal to our test function / the above ratio can be bounded from below using 
(|425|) and ([111 ) by 

N 3 - d 

C 2 QogN)*h ' (5 ' 29) 
where C 2 = Cin d h\ 2 . □ 

Clearly in dimension d ^ 3 the test function does not provide any information. 



5.2. Proof of Theorem |2.2| . Fix d ^ 2. By definition 

^+(0 2 log 2 ) 



When (ft is equal to our test function / the above ratio can be bounded from below using 
( |425|) and (p~2l] ) by 

(5.30) 



C 3 (logN)dh 



where C- 



a 



3 ClK d /TA 2 " 



□ 



6. Slow down of the Glauber dynamics in two dimensions 

In this section we will derive some consequences from the two dimensional upper bound 
on the inverse spectral gap for (3 > f3 c on the speed of relaxation of the Glauber dynamics 
to its equilibrium. In particular we will prove theorem ^L^. The notation will be that fixed 
in section 2 and section 3. 



is a 



6.1. A first finite volume bound. The first simple consequence of Theorem 2.1 
bound on the dynamical evolution of the test function ( 4.1 2| ) itself 

Proposition 6.1. For any N large enough, 

Vt > 0, ^ ((P t i+ ' N) ff) > »Uf) (-2*^^) I 1 " exp(-c A A0) , (6.31) 
where c\ is a positive constant depending on A. 

This result provides a first (admittedly weak) clue on the relaxation time of the dy- 
namics. Let us assume that the Markov process generated by is attractive (see |L| or 
[Maj ). This is the case if for example the flip rates were those of the Metropolis or of the 
Heat Bath dynamics. Let B~^{a) be the number blocks for which the spin configuration 
a in {iljWjv h as averaged magnetization smaller than — ^h- We set 

V N (a) =expf-^-0^(<7) 
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where K = (Mog AT) 1 / 7 and 6,7 are as in the previous section. 

Since Bj^ is a non increasing function of the spin variables, the monotonicity inequalities 
for attractive processes imply 

Pi + ' N, (* K )H > {^{(P^'^f)) 1 ' 2 > (^((P t (+ ' N) /) 2 )) 1/2 , 

where the symbol (— ) denotes the configuration in Wat for which all the spins are equal 
to —1. Inequality ( |6.31| ) implies that there is e > such that for all N large enough 

P^' M (-) >eJ[ p((2A- T -- £ )f- ( K). 

On the other hand, as in the derivation of ( [4 . 1 5|) , one can check that 

Therefore for any time smaller than ^ ^ N ^i K the quantity '^(Wat(o')) (— ) is much 
larger than the equilibrium expectation of *$>n : in the above special sense the system has 
not yet relaxed. 

Remark 6.1. It is important to observe that in the above reasoning we have never used the 
information that the logarithmic Sobolev constant is larger than ~ A^ 2 . Unfortunately we 
have not been able to establish anything like proposition i| for the entropy of [Pt + f) 2 
with the exponent t ( l ° s ^ replaced by t . 

Proof. We set 4> = f — fj^(f). The spectral decomposition of £^ implies 

, . poo 

/4 {{Pt' N U) ) = J <M*) exp(-2^) , 
where denotes the spectral measure associated to (j). By Jensen inequality, 

or *^)- 

By definition of the spectral measure 

f'OO f'OO 

li + N {4> 2 )= d»4,{0) and £%{f,f)= / 0<M<?). 
Jo Jo 

Therefore the bound on the spectral gap (see ( |5.29|) ) implies that for N large enough, 

Mjy {(?t' N) <P) 2 ) > ,W) exp (-2t^^) • (6.32) 
According to ( |4.15 ) , there is a constant c\ > such that 

/4(/) 2 ^(/ 2 )exp(-c A iV). 
The former inequality combined with ( |6.32| ) leads to 

Mjv {(Pt i+ ' N) ff) > ^(/ 2 ) exp (-2t^^j (1 - exp(-c A iV)) . (6.33) 
This concludes the proof. □ 
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6.2. Proof of Theorem 5T3. 



Proof. The first step is to reformulate the LHS of fl2.6|) in terms of mesoscopic variables. 
For any site x G if Z 2 we define Cr(i) to be the indicator function of the event that the 
magnetization in the box W>k(x) for the process ^(t) in Mk(x) is smaller than — Then 
we have 

A+(e(C WCo^))) < A + (lE(3,G B K (0), aV(t)^5f(t) 

< K 2 We(^)^oW 



where we used the invariance by spatial translation in the last inequality. 

Let N be a large integer, choose as usual the mesoscopic scale K = (MogiV) 1 / 7 and let 

L = jf. By repeating the previous computation on a coarse grained level, we get 

a+(e( £ cw/ £ crw)) < A + (E(3f G w iV nKz 2 , cm/crw 



lV (fi(^(t) + ft®)) . 



ieWivn^z 2 



Let now 



The previous results imply 

Aw (EK(t) / ogf(t))) > (E(B^(o-?) + B N (at))) . (6.34) 

In the second step, we are going to decouple the estimates of the joint process. The main 



physical idea was already contained in the Fisher, Huse paper [HF] and it goes as follows. 
We force one large droplet of the minus phase of radius ~ N, around the origin in e.g. 
the initial distribution of crj 7 , by paying a price ~ exp(— r* N). This droplet should relax 
only in a time scale proportional to its initial area and therefore, if N = A\ft with A large 
enough, the distribution of B^(a^) at time t given the above initial unlikely event should 
be quite different from that of Bjj(af). Apparently in order to carry rigorously the above 
program one needs a much more precise control on the life time of a droplet than what 
we have been able to obtain. Actually that is not true and all what we need is something 
not more precise than proposition [O] (see Lemma |6.1| below) . 

From a technical point of view it is convenient to force the droplet of the minus phase 
inside fJ- + (rj) in a "soft" way by simply inserting our test function / 2 defined in ( 4.12Q with 
N^Vt. 
We write 

f,+ (E(B N (a?)^B N W))) > 
> ^+(/ 2 )exp(-2AiV) J d^ 2 ( V )dv+(u)E{B N (a?)^B N (a?)) > 

> exp(-2r*A0 j df'? (rj) d^{u) E(B N (a r t ') + B~ N {a?)) (6.35) 

Let a be a parameter in (0, 1) which will be fixed later on. Then 

&{B- N {a*) + B N (a?)) > E{B N (^) > aL 2 ; B N (a?) ^ aL 2 ) 

> E(0-(a?) > aL 2 ) - E(B^(a?) > aL 2 ) , 
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where E refers to the marginal of E, i.e. to the usual Glauber dynamics. Since the measure 
fi + is invariant with respect to the Glauber dynamics, we can write 

fj+ {E(B N (^) > aL 2 )) = ii + (B N {o) > aL 2 ) ^ exp(-C a N) , 

where the final estimate follows from the theory of equilibrium phase coexistence (see 
Propositions |3T3| and |3.1| ). In conclusion 

It is at this stage that we are going to use the information on the spectral gap. The 
necessary dynamical estimate is provided by the following Lemma which will be derived 
later. 

Lemma 6.1. We fix a such that the parameter (2A(a — 1) + r*) is negative. Then, for 
N large enough, the following inequality holds 

L (logJV)^ 



Vf 



> 0, J d^f\n) E{B N (a?) > aL 2 ) > \ exp(-t - exp(-c a , A iV) . 



N 



where c Q; a > 0. We recall that L = 
The previous Lemma implies that 



exp(-2r*iv) Q e xp(-t^^) - exp(-c a , A iV) - exp{-C a N) j . 

By choosing N = \ft (logt) K , we finally derive for t large enough 

(i+(E(B^)^B^))) > e W (-2r*Vi(logtr-3Vt). 



□ 



Proof of Lemma 6.1. The proof relies on the dynamical estimate of proposition 6.1 
Let 



*t = W(p t /) 



E(/(a?)l 



)+E(/(^)l B - (CT?) > ai2 ) • (6.36) 



Thus from the estimate ( fhl4|) and the FKG inequality we see that for e small enough, a 
such that (2A(a - 1) + r*) = 5 a < and N large 

E (/K") ^ (fl? )<^) < exp(aAiV) = exp ((2X - r* - * a )y) < 

< exp(-^iV)^+(/2) < exp( _^ iV )7^T(72y. 

Plugging the above inequality in fl6.36j ), we get 



E (/(°7) +^(/ 2 )ex P (-f n), 



(6.37) 



with 5 a > 0. 

In the pure phase // + , the estimates obtained in subsection 4.1 for the variance and the 
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Dirichlet form of / hold (see Proposition |6.2| below). Thus Proposition 6.1 is also valid for 
an unbounded region and for N large enough, we get 

Vt > 0, % > /z+(/ 2 ) exp (-2t {l ° S * r ^ (1 - exp(-c A iV)) . 



Combining the previous inequality with ( 6.371 ), we get by using Cauchy Schwartz in- 
equality 

M + (E(/ 2 K))E(l s _ K)>aL2 )) > M + (/ 2 )(iexp(-2ti^l!)-ex P (-c a iV)) . 
The reversibility of the dynamics ensures that 

^ (e(/ 2 (^))e(i b - K) ^ l2 )) =m + (/ 2 (i)E(i B - (4) ^ i2 )) • 

This concludes the Lemma. □ 
Proposition 6.2. In dimension d = 2, for any (3 > (3 C then 

ViV, M +(/ 2 ) - M + (/ ) 2 > C p^M + (|V/| a ) , 



where f is the test function introduced in ( 4.12| ) . 



Proof. The upper bound ( 4.25| ) on the Dirichlet form is unchanged under [i + since it 



boils down to estimating the number of bad blocks in the region Wtv by using Bernoulli 
percolation. The lower bound fl4.14j ) holds also for /u + (/ 2 ) because it involves only the 
computation of an event localized in W^v- 
Thus it remains only to check that 

limsup — : log fi + {f) ^ 0. (6.38) 

iV->oo N 

Let D = [— R, R] 2 C M. 2 , where R will be chosen large enough; in particular such 

that W 2 C [— R/2, R/2] 2 . Let Cn be the set of spin configurations which contain a *— 
connected circuit of + spins inside [— N, N] 2 \ [— N/2, N/2] 2 separating ([— N, N] 2 ) c from 
[- N/2, N/2} 2 . As (3 > (3 C and d = 2, there is cp > such that 

fi + (C c N ) ^ exp(-cpN) . 

By choosing R such that Rep > 2A, we get 

V + (f) ^ V + (f lc NR ) + V + (C C NR ) exp(AiV) < fx+(f 1 Cnr ) + exp((A - c p R)N) . 

Conditionning with respect to the + circuit which is the closest to ([— RN, RN] d ) c and 
then using the fact that / is non-increasing, we obtain by FKG 

where /ji^r denotes now the Gibbs measure in {—NR, . . . , NR} 2 with + boundary condi- 
tions. 

At this point, we proceed as in subsection 4.1. The only difference is that the estimates 
are in iJ-(D) instead of L 1 (W 2 ). This implies 

limsup -j- log fj,+ R {f) < sup { -Wfs(v) + X\{v = -l}nW 2 \}. 



ueBV(D,{±l}) 
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To derive ( 6.38| ), it remains to check that the RHS is negative. Either |u| ^ 1 in which 



case, we get 

-VW^(u) +\\{v = -1} nW 2 | < -T*\v\ {d ~ 1)/d + A|W 2 | s: -r* + A<0. 



or \v\ < 1 and (]4. 19 ) applies 

-Wp(v) +\\{v = -i}nw 2 | < -t*\v\ +X\v\ < 0. 



□ 



Remark 6.2. As a consequence of the proof of Proposition \6. ij , we see that the inverse 
of the spectral gap associated to the Glauber dynamics in the cube [— N, N] 2 grows faster 



than x 



(logJV)«- 

7. A ONE DIMENSIONAL BIRTH AND DEATH PROCESS FOR THE DROPLET EVOLUTION 

In this section we discuss a simple one dimensional toy model which mimics the random 
evolution of the volume of a droplet of the minus phase in a large cube of side N in Z d 
under a Glauber dynamics with plus boundary condition. The model goes as follows. 
Let a := and consider a birth and death process on the integers A := {0, . . . , N d }, 
reversible with respect to the measure 

MO) : = 7? exp(-x a ) 

and with birth rate b(x) = (x V l) a , x < N. By reversibility the death rate d(x) is given 
by 

d(x + 1) := (x V l) a exp((x + l) a - x a ), x>0 

One easily checks that the drift given by b(x) — d(x) is negative and proportional to ai 2a_1 
for large x. The connection with the evolution of a large droplet of the minus phase under 
the Glauber dynamics with plus boundary condition in a large cube of side N in 7L d is as 
follows. The variable x represents the volume of the droplet at time t which is assumed 
to form a compact set without holes. The rate b(x) should then be interpreted as the rate 
with which a plus spin just outside the boundary of the droplet flips to minus one and gets 
attached to the droplet while the rate d{x) represents just the opposite process in which a 
minus spin at the boundary flips to plus one and gets detached from the droplet. Clearly 
both these rates should be proportional to the size of the boundary which, for roundish 
shape, is of order of x a . Finally the drift comes from the reversibility condition together 
with the fact that the equilibrium distribution of the droplet volume should behave like 
the measure fi(x) above according to the results of section ||[ Quite nicely the drift one 
gets out of these natural hypotheses is of the same order of that prescribed by an evolution 
by mean curvature 

± x i/d = L 



dt x 1 / d 

Our goal now is to compute the precise asymptotic as N — > oo of the inverse spectral gap 
S(N, d) and logarithmic Sobolev constant L(N, d) of the above process in order to test 
the accuracy of the bounds proved in section |2.3|. 



Theorem 7.1. For any d > 2 there exists a positive constant k = k{d) such that 
(i) ^ ^ L(N,d) N 2 , yd ^ 2 



18 



T. BODINEAU AND F. MARTINELLI 



(ii) f < S{N, 2) < kN 

(iii) S(N,d) < k, Vd ^ 3 

Proof. We apply the method of Hardy inequalities envisaged in |Mi| |) in order to compute 
sharp upper and lower bounds on the quantities of interest. We begin with the inverse 
spectral gap and define 

B+(i): ^(i™)|> w 

B:=inf(B + (i)VS_(i) 

The measure /i is extended on Z by setting n(x) = if x {0, . . . , iV d }. Then we have 
(see Proposition 1.3 of (||Mi||) 

— «S 5(iV, d) < 4£ 

Part (ii) and (iii) of the theorem follow at once from the simple estimates 

^(y)^ 1 -" exp(-x Q ) 

y>x 

Y 7 !,/ \ ~ ^- 2a exp(x") (7.39) 

where A ~ means that there exists a universal constant such that f ^ 4 ^ We 
get in fact that for i E {0, . . . , iV d }; B + {i) ^ N ior d = 2 and B + (i) ^ uniformly in iV 
for ci ^ 3, while B-(i) ~ i 1_2a exp(i a ) for any <i. Notice that B + {i) = oo if i < and 
= oo if i > N d . 

We now turn to the analysis of the logarithmic Sobolev constant. We define 
A+(i) :=sup( V - 1 ) log( = — — 7-r)yZKy) 

i-1 . j 

A_(») := sup( V / \ ) lQ g(^ rr) y2Kv) 

x<i^^ x Ky)Ky) J K E y < x Ky) J ^ 



A:= inf(A + «VA_« 
Then we have (see Proposition 3.1 of [[M]]) 



—A < L(N,d) < 2(L4 
20 v 7 



and part (i) follows at once from the bounds ( ^.39 ). 
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